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Abstract 

Using molecular dynamics simulations we investigate the translational dynamics of particles with 
dipolar interactions in homogenous external fields. For a broad range of concentrations, we find 
that the anisotropic, yet normal diffusive behavior characterizing weakly coupled systems becomes 
anomalous both parallel and perpendicular to the field at sufficiently high dipolar coupling and field 
strength. After the ballistic regime, chain formation first yields cagedike motion in all directions, 
followed by transient, mixed diffusive-superdiffusive behavior resulting from cooperative motion of 
the chains. The enhanced dynamics disappears only at higher densities close to crystallization. 

PACS numbers: ??? 



There is a strong current interest in colloidal soft-matter systems displaying anomalous 
dynamic behavior beyond simple diffusion. A paradigm example are colloidal hard spheres 
which exhibit a (repulsive) glass transition in near-perfect agreement with Mode Coupling 
theory [l]. Recent attention focuses, on one hand, on the anomalous slow dynamics in a 
)road variety of systems such as low-density colloidal gels Q, cells Q, driven granular fluids 
4], and disordered solids j^j. In many cases the reduced mobility is accompanied by dynamic 
heterogeneities and cooperative motion as recently shown, e.g., in colloids js, 7] and in liquid 
crystals 8|. On the other hand, there is increasing interest in systems exhibiting enhanced 



mobility such as rods between obstacles Jq] and activated systems such as catalytic nanorods 



10( or self-propelled bacteria [11]. 

In this letter we discuss the translational dynamics of dipolar systems which, as we will 
demonstrate, can exhibit both, slow and enhanced mobility. Prominent examples of dipole- 
coupled systems are suspensions of magnetic nanoparticles (ferrocolloids), polarizable col- 
loids [12[ or colloidal "molecules" (e.g. dumbells) with partial charges [131 ]. but also elec- 
tric field-driven living cells jl^j]. The strong current attention received by dipolar systems 
(or, more generally, systems with directional interactions) is born by their importance for 



sel 
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-assembly processes towards novel classes of equilibrium and non-equilibrium materials 



161 ] . In zero field, the anisotropic interactions between permanent dipoles can induce 
network formation |l7j including the possibility of an equilibrium gel 18] . On the contrary, 
the dynamics in an external field appears far less understood, although this situation is 
clearly most relevant for many practical applications of dipolar systems 19[. Previous ex- 
perimental and theoretical studies indicate that the single-particle motion under field can be 
strongly anisotropic 20[ and even sub-diffusive (in lateral directions) due to " caging" of the 
particles in the chains formed at sufficiently hi gh c oupling strength 2ll423j|. However, only 
recently dynamic light-scattering experiments 24J reveal, in addition, strong evidence for 
super diffusive motion and dynamic heterogeneity. Motivated by these exciting and seem- 
ingly somewhat contradictory findings, we have performed extensive molecular dynamics 
(MD) simulations of a simple dipolar model system. Analyzing separately the microscopic 
motion parallel and perpendicular to the field for a broad range of parameters we find, for 
the first time, a coupling strength-induced crossover from normal to anomalous dynamics. 
The latter is characterized by a step-wise behavior involving sub- diffusion, directed motion 
and dynamic heterogeneities. 
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FIG. 1: (Color online) a) MSDs for A = 3 at H* = and 10; b) simulation snapshot at A = 3 and 
H* = 10; c) anisotropy ratio as function of dipolar coupling strength. 

Our model fluid consists of dipolar soft spheres (DSS) of diameter a with embedded, 
permanent point dipole moments Hi. The total pair interaction at distance r = |r| consists 
of a (truncated and shifted) soft-sphere repulsion, «ss( r ) = 4e(cr/r) 12 , and the long-range, 



anisotropic dipole-dipole potential, u DD (12) 



[Hi ■ n 2 ~ 3 (jii • r) (n 2 ■ r) /r 2 ]. We em- 



ploy constant-temperature MD simulations with N = 500 — 1372 particles. Thus, we focus 
on the direct (conservative) effects on the dynamics rather than on those stemming from a 
solvent, in accordance with other recent studies (see, e.g., 18]). The dipolar interactions are 
treated using the Ewald method with conducting boundaries (see 25J for further details). 
Simulations are carried out at a reduced density p* = per 3 = 0.05 (which is typical for the 



experiments in |24|) and a reduced temperature T* = k-^T/e = 1.35 (with ks and T being 
Boltzmann's constant and true temperature, respectively). The dipole moment is chosen 
such that the resulting coupling parameter A = fi 2 /kBTa 3 has values 1 — 7, consistent with 



the range considered in recent ferrofluid experiments [2j, [26|. The fluids are subject to a 
homogeneous magnetic field H = He z with fields strengths up to the (still realistic 261 ]) 
value H* = fiH/k B T w 100. 

One key quantity of our analysis is the translational mean-squared displacement (MSD), 
Ar 2 (t) = ((r(£) — r(0)) 2 ). Representative MSDs for a moderate coupling parameter (A = 3) 
are shown in Fig. (T^i). For H > we differ between the MSDs parallel and perpendicular 
to the field direction. This is sensible due to the structural anisotropy resulting from the 
formation of elongated clusters along the field, see Fig. [lb). Despite these (small) clusters, 
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FIG. 2: (Color online) a)-b) MSDs for A = 7 and H* = 0, 10, 100. b) includes data for two system 
sizes; the box indicates the sub-diffusive regime. c)-e) corresponding snapshots. 

all of systems have "normal" translational dynamics in the sense that there is an initial 
ballistic regime, Ar|^_g (t) oc t 2 , after which the MSD directly goes over into a diffusive regime 
Arjj/jg(t) oc t. The latter, Einstein-like behavior allows us to define diffusion constants Dq 
(for H = 0) and D», D± (H > 0) in the standard way. Of particular interest is the anisotropy 
ratio A = [D\\ — D±) /Dq, results for which are plotted in Fig.dfc). At the coupling strengths 
considered, motion perpendicular to H is "faster" than along the field, i.e., A < (consistent 



with earlier simulation results [201]). despite the fact that the clusters form in z-direction. 
Thus the lateral translational fluctuations of the particles are still so large that their diffusive 
behavior resembles disks than rods (where one expects A > 0). However, as Fig. [T^) clearly 
indicates, there is a crossover towards different dynamic behavior as A increases beyond ~ 5. 

Typical MSD's for this strongly coupled regime are plotted in Fig. [2] along with corre- 
sponding simulation snapshots (A = 7). At zero field, the particles arrange into a globally 



isotropic network of worm- 
ing dipolar systems 



17 



ike chains [Fig. Efc)], as it is typical for dilute, strongly interact- 
18l. The corresponding MSD still behaves "normal" in the sense 

n 

described above, consistent with MD results for dipolar dumbbells [18j at comparable values 
of A. With field, long chains form which become more and more straight upon increase of H. 
We first consider the case H* = 100 [see Figs. [2b) and e)]. After the initial ballistic regime, 
the motion of a particle significantly slows down due to the confinement by its nearest neigh- 
bors, yielding a sub-diffusive regime in the MSD characterized by a time dependence t 7 with 
7 < 1. This holds not only in z-direction (where sub-diffusion occurs for 0.6 < t* < 1 with 
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FIG. 3: (Color online) a), b): Distinct- van Hove functions (A = 7, H* = 100). c), d): Trajectories 
of two initially neighboring particles (insets show a larger time scale). 

711 ~ 0.61), but also perpendicular to H (0.3 < t* < 3.3, 7^ ~ 0.69), indicating pronounced 
spatial-temporal correlations in all spatial directions. We characterize these correlations by 
the distinct parts of the van Hove correlation functions G\(z,t) and G^(R,t) with z and R 
being the longitudinal and lateral distance, respectively. Results are plotted in Figs. 13^,) and 
b). For t — > the functions G|| coincide with the corresponding static pair distribution 
functions. The data for G\ at t > clearly reflect long-range correlations, and thus, a 
"trapping" of particles along the chain, which starts to die out only at times beyond the sub- 
diffusive regime of the MSD (i.e., t* > 1.1). The exponent 711 ~ 0.61 within the sub-diffusive 
regime is close to the value of 0.63 found in the MSD of interior beads in linear polymer 
chains {2^ prior to crossover to free diffusion. In perpendicular directions [see Fig. [3b)], the 
correlation "hole" at t — > reflects that positions directly besides a particle are disfavored 
due to the repulsive dipolar interactions associated to this configuration. At later times, 
however, we observe development of a peak at R — > 0, indicating that the original position 
of a particle is replaced by another particle. These strong and long-living correlations ex- 
plain the sub-diffusive behavior of the perpendicular MSD, characterized by an exponent 
7x = 0.69. A similar exponent (7 w 0.75) has been found in experiments probing lat- 



eral diffusion in magnetorheological systems 2l|, |22| involving chains of induced (and thus 



perfectly aligned) dipoles. This similarity already indicates that the strong suppression of 
orientational fluctuations, i.e., a strong field, is an important ingredient for the observed 
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sub-diffusive behavior in our present system consisting of permanent dipoles. Indeed, the 
MSD's at the somewhat weaker field H* = 10 [see Fig. Hk),d)] do not reveal a sub-diffusive 
regime. Instead, one merely observes a slowing down (7y = 1.13, -yj_ = 1.12) as compared to 
the initial ballistic regime. 

However, both anisotropic systems (H* = 100 and H* = 10) display, at A = 7, a further 
intriguing feature not present in weakly coupled systems considered in Fig. [TJ Specifically, 
we observe in Figs. [2k)-b) a "kink" in the MSDs, which separate the intermediate time range 
discussed above and a later time regime where the MSDs grow faster than linearly in time. 
This behavior strongly contrasts that in other disordered systems (see, e.g., j^]), where the 
sub-diffusive regime, if present, changes directly into true diffusion (7 = 1). In our present 
system, the initial exponent within the "fast" regime (after the kink) lies between 1.5 and 
1.8; however, these values decrease in time and also strongly depend on H*, as does the time 
related to the kink. We conclude that there is a (transient) superposition of diffusive and 
ballistic (7 = 2) behavior, with a tendency towards normal diffusion (visible particularly in 
the lateral part) at large times. In fact, at the low density considered (p* = 0.05) we did not 
reach the true diffusive regime particularly for the longitudinal part in the total simulation 
time (contrary to what is found at larger p, see Fig. [5]). To test possible finite-size effects we 
have repeated the simulations with N = 1372 particles, the changes [see Fig. [2b)], however, 
being only marginal. 

A direct illustration of the superposition of diffusive and ballistic motion is provided 
by typical trajectories plotted in the bottom of Fig. [31 In ^-direction [3b)], one clearly 
recognizes small random- walk- like movements of the particles around a mean path consisting 
of essentially straight stretches, the latter giving rise to the "fast" regime of the parallel 
MSD [see Fig. [2]. Similar random movement around a mean path occurs perpendicular to 
H [Fig. [3H)], where, however, the stretches are shorter and, consequently, the growth of the 
MSD (after the kink) is somewhat less pronounced. Moreover, in both directions, the motion 
of two particles which are nearest neighbors at t = is correlated over an extremely long time, 
consistent with the behavior of G\ (and ) shown in the top of Fig. [3j These observations 
suggest that, within the fast regime, the particles move collectively as chains (indeed, at 
the large field considered, there are essentially no unbounded particles). Moreover, the 
stretches in the trajectories indicate that the entire chains can move quasi-independently 
(i.e., ballistically) over a finite time interval before they are slowed down due to interactions 
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FIG. 4: (Color online) a), b) The functions Gs at t* = 20 with corresponding instantaneous 
data. The insets show the non-Gaussian parameters a\\f±\(t). c), d) Snapshots with particles 
colored according to the code in a) and b). 

with particles in other chains. 

The occurrence of cooperative motion of the particles as chains is further confirmed by 
the behavior of the self parts of the van Hove functions, Gs (z, t) and G^(R, t) measuring the 
probability that a particle moves over a certain distance in a time t. Exemplary results for 
time t* = 20 (i.e., within the "fast" regime) are shown in Fig.HJ Particularly the parallel part 
strongly deviates from the usual Gaussian shape characterizing diffusive systems. Similar 
deviations occur at other times (and in the perpendicular part), yielding large values of 
the "non- Gaussian" -parameter a(t) = ((Ar) 4 )/((Ar) 2 ) 2 — 1 (see insets in Fig. H]). At t* = 
20, specifically, the buckled structure of Gl [Fig. @Ji)] suggests that the system consists of 
several portions ("populations") of particles. To better identify these populations we have 
resolved the instantaneous data for G s with respect to the actual particles in the system 
as illustrated in Fig.Hfc),d). One clearly recognizes that particles related to one "population" 
in Gs belong to the same chain in the actual configuration. Also, comparison of Figs. Hfc) 
and d) shows that chains which are fast in ^-direction move more slowly in perpendicular 
directions and vice versa. 

Taken altogether, we interpret the "kink" in the MSDs (and thus, the superdiffusive 
behavior), as a consequence of the motion of the entire chains, which over certain times, 
interact only weakly with particles of other chains. This interpretation prompts the ques- 
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FIG. 5: (Color online) MSDs a) parallel b) perpendicular to H at various densities p* (A = 7, 
H* = 100, N = 864). Part c) is a snapshot at p* = 0.8 



tion what would happen at higher densities, where chain-chain interactions should be more 
prominent. To this end, we plot in Fig. |5] results for MSD's in the range 0.05 < p* < 0.8. 
The simulations have been performed at the same strong coupling conditions as before, 
yielding essentially percolated chains in z-direction. An inspection of the data in Fig. re- 
veals that the "kink" indicating onset of (transient) super diffusive motion along H survives 
up to p* ~ 0.3 in field direction; at higher p* one observes normal dynamics where the ini- 
tial ballistic regime directly changes into diffusion. Perpendicular to H, the disappearance 
of the kink occurs already p* ~ 0.1; at p* > 0.3 plateau-like regions in the MSD appear. 



Indeed, at high densities one expects 28J formation of body-centered-tetragonal structures 
with neighboring chains being out of registry, a tendency which is visible in the snapshot 
in Fig. |5fc). Clearly, this lateral ordering induces additional dynamic effects not present at 
small p*. 

To summarize, our MD simulation results demonstrate that chain formation of strongly 
coupled, dilute or moderately dense dipolar fluids in external fields gives rise to anomalous 
translational dynamics. The latter involves both, sub-diffusive behavior at intermediate times 
and mixed diffusive-superdiffusive behavior at later times, consistent with recent ferroffuid 



experiments 



24J . According to our analysis in the superdiffusive regime, the particles move 



collectively as chains, a behavior which appears as dynamic heterogeneity on the single- 
particle level. From a fundamental point of view, a particularly intriguing result of our 
study is that the systems may display enhanced dynamics already in equilibrium, cont rary 
to active (bio-) materials such as self-propelled bacteria and catalytic nanorods 10) . 
which are intrinsically out of equilibrium. From an application perspective, we expect our 
results to be relevant not only for ferrofiuids, but for the much broader class of colloidal 
systems where external fields induce directed self-assembly, examples being patchy colloids 
such as metallodielectric (Janus) particles 15[ and living cells suspensions 14|. Clearly more 
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detailed experiments are needed to confirm our findings, such as the "kink" separating the 
dynamic regimes in the MSD. On the theoretical side, points deserving further attention are 
the role of orientational fluctuations at moderate fields, as well as solvent (hydrodynamic) 
effects which are so far neglected. Work in these directions is in progress. 

Financial support from the German Science Foundation via grant KL1215/6 is gratefully 
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